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Quantum Monte Carlo (QMC) 


• Quantum Monte Carlo provides reliable solutions to quantum 
many-body problems. 

• It can also be used in classical optimization problems as an alternative 
to simulated annealing. 

• Whether quantum tunneling can be simulated efficiently by quantum 
Monte Carlo? 


< □ ► « S? ► « 1; ► < I ► ^ 


Zhang Jiang (QuAIL) 


Instantons in quantum annealing 


December 16, 2015 


3/22 





Why QMC? 


• Recent Google result (arXiv:1512. 02206) showed that there is no 
asymptotic speed up when the D-Wave quantum annealer is 
compared to QMC, although a constant factor 10 8 is observed. 

• There exists an analogy between the tunneling decay of quantum 
systems and classical escape-over-a-barrier problem. 

J.S. Langer, Theory of Condensation point, Annals of Physics 41, 108 (1967). 

M.Buttiker and R.Landauer, Nucleation theory of overdamped soliton motion, 
Phys. Rev. Lett. 43, 1457 (1979). 
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What is QMC? 


Path-integral Monte Carlo works by Trotterizing the partition function in 
imaginary time 






Tunneling in spin systems 


We consider the mean-field model where the system Hamiltonian is 
symmetric with respect to permutation of individual spins 

H = —NVm x — Ng(m z ) 

1 N 

19)01 = /v a< ? ’ oi = x,y,z 

v 1=1 

Here g(m) is a nonlinear term that allows for co-existing local and global 
minima for m E (-1,1). 

The WKB Hamiltonian is 

Hwkb ( m, p) = —21 " N \J £ 2 - m 2 cos p — Ng(m ) , 

where £ = 2S//V E (0, 1) and m = m z E (—£,£). 
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QMC probability functional 


We define the state vector 


H r ) = {<71 (r), . . . , <7/v(r)}, ct(0) = a(fl) 


The probability of a state vector is 


P[ct(t)] = Z 1 exp [—(3NE[(j(t)]\ , 
f [7(r)] = -- g[m{T)]dr - -y- 


«M r )] • 


The function k equals to the number of times ctj(t) changes its sign. 
Order parameter: the total magnetization 





QMC probability functional in reduced space 


There is a Gibbs probability measure P[m(r)] = Z 1 e N i 3F l m ( T )] for the 
magnetization order parameter m(r) ( Bapst , Semerjian, 2012) 

F [ m ( T )] = ^ j o l m ( T )g'( m ( T )) ~ g( m ( T ))] dr ~ Jj log A[g / (m(r))] 

Here, the functional A[A(r)] is 

A[A(r)] = TrA' / 3 ’°[B(r)], K T2 ’ T1 = T + e~ f n drH ° {r) 

H 0 {t) = -B{t) ■ a, B(t ) = (T, 0, A(r )) , 

where a = (a x , a y , <j z ) is vector of Pauli matrices. The propagator K 
corresponds to a spin-1/2 evolving under the magnetic field B(t). 
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Kramers escape problem 





The system reaches the transition state via thermal fluctuation. Then with 
probability ~ 1/2 it moves toward the global minimum. 

« fit ► «!;► < 1 ► ^ 


'O 0,0 


Zhang Jiang (QuAIL) 


Instantons in quantum annealing 


December 16, 2015 


9/22 





Kramers escape in QMC 


QMC samples paths a(r, t ). When the path m(r, t) moves from the local 
minimum to the global minimum by fluctuation, it has to go through the 

transition state m z (r) 

Wqmc OC e~ l3NAF , A F = F[m z (r )] - F(m 0 ) . 


Here m z (r) is the saddle point of the functional F that satisfies the 
equation 


SF i m ( T )] 

Si7i(t) 



m z { 0) 


m z{(3) ■ 
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Instantons in QMC 


Variational equations SF = 0 take the form 

6 log A(A(r)) dU[m z (r)\ 

m ' (T) = JA(r) ’ A(T) = </m 2 ' 

We introduce vector of magnetization components 

, . Tr[KP’ T aK T '°] 
m(r) = TrK8l , 

Optimal trajectory is a classical rotator in nonlinear potential 
dm(r) _ _ 2 .dHo[m(T)\ 


dr 


3m 


x m(r) 


'H 0 [m\ = -r m x (r) - L/[m z (r)] 
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Instantons in QMC II 


Two integrals of motion 

Tfo[m] = e, m(r) • m(r) = 
Then the solution can be written in the following form: 



The equation for m z (r) is identical to that of the WKB instanton 
trajectory. 
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QMC for Curie Weiss model 







Optimal value of total spin (f) Logarithm of escape rate In) 


Comparison of QMC and 
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Future work 


• Quantum Monte Carlo with open boundary condition 

• Calculate the prefactor for QMC and quantum annealing 

• More general spin coupling 
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